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Abstract 

We combine a recent mapping of the Anderson-Mott metal-insulator transi- 
tion on a random-field problem with scaling concepts for random-field mag- 
nets to argue that disordered electrons near an Anderson-Mott transition 
show glass-like behavior. We first discuss attempts to interpret experimental 
results in terms of a conventional scaling picture, and argue that some of the 
difficulties encountered point towards a glassy nature of the electrons. We 
then develop a general scaling theory for a quantum glass, and discuss criti- 
cal properties of both thermodynamic and transport variables in terms of it. 
Our most important conclusions are that for a correct interpretation of ex- 
periments one must distinguish between self-averaging and non-self averaging 
observables, and that dynamical or temperature scaling is not of power-law 
type but rather activated, i.e. given by a generalized Vogel-Fulcher law. Re- 
cent mutually contradicting experimental results on Si:P are discussed in the 
light of this, and new experiments are proposed to test the predictions of our 
quantum glass scaling theory. 

PACS numbers: 71.30.+h, 05.30.-d, 75.10.Nr 
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I. INTRODUCTION 



The metal-insulator transition that is observed in doped semiconductors and other dis- 
ordered solids is not fully understood, despite almost twenty years of intense experimental 
and theoretical effort s.0 There is strong evidence for both disorder and electron-electron 
interactions to play an important role at this transition, which is called an Anderson-Mott 
transition (AMT) to distinguish it from a disorder dominated pure localization or Ander- 
son transition on one hand, and from a correlation dominated pure Mott transition on the 
other hand.! This interplay between disorder and interactions makes the AMT a very hard 
problem in Statistical Mechanics. 

On the theoretical side, until very recently virtually all approaches studied the problem 
in the vicinity of two dimensions (d=2) by generalizing Wegner's theory for the Anderson 
transition.! These theories^ have led to a classification of the AMT into various universality 
classes that depend, inter alia, on the presence or absence of spin-orbit scattering, magnetic 
impurities, magnetic fields, etc. For most of these universality classes, perturbative renor- 
malization group methods lead to a critical fixed point in d = 2 + e dimensions, and standard 
critical behavior with power-law scaling is found. However, the framework of these theories 
does not allow for an order parameter (OP) description of the AMT, and does not lead to 
a simple Landau or mean-field theory. As a result, the physics driving the AMT remains 
relatively obscure in this approach, compared to standard theories for other phase transi- 
tions. An alternative line of attack has recently been explored by the present authorsJUl. In 
these papers we showed that the same modelQ that was used for the 2 + e expansion allows 
for an OP description of the AMT with the tunneling density of states (DOS) as the OP, 
and for a simple Landau theory that yields the critical behavior exactly above the upper 
critical dimension d^ = 6. Furthermore, it was shown that the problem has random-field as- 
pects and is closely related to a random-field Ising model. The structure of that theory also 
suggests that for certain parameter values, in particular for weak effective electron-electron 
interactions, the OP driven AMT can be pre-empted by a different kind of metal-insulator 
transition. The most obvious candidate is the Anderson transition, where the DOS or OP is 
uncritical. These results lead to the suspicion that in large parts of parameter space impor- 
tant physical features of the AMT have been missed both in the low-dimensional theories, 
and in the interpretations of experiments that were based on these theories. 

On the experimental side, the best studied systems are doped semiconductors, and some 
of the most detailed and careful experiments have been done on Si: P. In what became a 
benchmark experiment in the field, Paalanen, Rosenbaum, Thomas and collaboratorslHill (to 
be referred to as the Bell experiment) performed transport measurements at temperatures, 
T, down to 2.7 mK, and used a stress tuning technique to drive a barely insulating sample 
through the metal-insulator transition. These experiments concluded that the critical P 
concentration, n c , in this system is close to 3.7 x 10 18 cm -3 , and that the conductivity, 
extrapolated to T = 0, vanishes with a critical exponent s ~ 0.5, with error bars of about 3% 
for n c and about 14% for s. The measured value of s proved hard to understand theoretically, 
and the critical behavior of Si:P is still considered enigmatic.0 More disturbingly, however, 
a similar experiment on the same systemS (to be referred to as the Karlsruhe experiment) 
has recently produced results that are inconsistent with those of Refs. |3|-|I~0|, viz. n c close to 
3.5 x 10 18 cm -3 , and s ~ 1.3. These disagreements are far greater than the error bars quoted, 
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and have not been settled between the respective experimental groups.El Measurements 
of the Hall coefficient have also led to mutually contradicting results. Dai, Zhang, and 
SarachikEl have reported evidence for a divergent Hall coeficient in Si:P, in contrast to 
previous results^! that had found the Hall coefficient to remain finite. Finally, different 
experimental results have been obtained for the crossover exponent of the (longitudinal) 
conductivity in an external magnetic field,0 but in this case no direct comparisons on the 
same system have been done. 

The fact that the most careful experiments on the best studied systems lead to contradic- 
tory results is, taken at face value, extremely discouraging. An important question is what 
the source of these discrepancies is. One relevant consideration is the temperature range cov- 
ered by the respective experiments, and the lowest temperature reached. This is important 
since any measurement of static critical exponents at the AMT involves an extrapolation to 
T = 0. For instance, the discrepancy concerning the Hall effect results has been blamed in 
Ref. H on the earlier experiments not having reached sufficiently low temperatures. In the 
case of the conductivity exponent in Si:P, however, extrapolation problems by themselves 
are not sufficient to explain the disagreement. References [|,|9] and [L^ agree that the ex- 
trapolated conductivity shows a 'tail' at phosphorus concentrations n < 3.7 x 10 18 cm~ 3 , 
but disagree about whether or not this tail contains the salient physics. The Bell group 
observed strong sample-to-sample variations of the conductivity in the tail region at the 
lowest temperatures, and concluded that the tail was due to sample inhomogeneities and 
should be discarded. The Karlsruhe group, on the other hand, claims that the tail repre- 
sents the asymptotic critical region of the AMT. The fact that no strong sample-to-sample 
fluctuations were observed in this case should, however, not be overweighted: The lowest 
temperature reached in the Karlsruhe experiment was T = 62.5 mK, while the strong fluctu- 
ations in the Bell data set in only at lower temperatures. Also, sample-to-sample variations 
in the T-dependence of the conductivity does seem to set in just at the lowest T reached 
in the Karlsruhe experiment, a phenomenon attributed in Ref. [T2| to 'thermal decoupling', 
i.e. problems in reaching and maintaining equilibrium. To summarize the experimental 
situation regarding the conductivity in Si:P, one might say that unusual features, variably 
described as 'rounding', 'smearing', 'thermal decoupling', etc. were observed at low temper- 
atures close to the critical point. These anomalies became stronger upon approaching the 
critical point and lowering the temperature, and the main discrepancy between the Bell and 
Karlsruhe experiments can be traced to different assumptions concerning their significance, 
i.e. whether or not the 'tail' should be taken seriously. 

Apart from the conductivity, unusual behavior has also been observed in thermodynamic 
properties of doped semiconductors. Both the magnetic susceptibility, Xm, and the specific 
heat, cy, show a pronounced non-Fermi liqid behavior JllHl! This behavior is observed near 
the transition as well as far away from it, is not obviously related to anv critical phenom- 
ena near the AMT, and is usually explained in terms of local moments.!^ However, Bhatt 
and Fisher§ have argued that once local-moment/local- moment interactions are taken into 
account, one obtains singularities in \ m and cy that are significantly weaker than those 
observed experimentally. This suggests that phenomena other than local moments might 
contribute to the observed anomalies in the thermodynamic susceptibilities. 

If one insists on a conventional theoretical interpretation of the conductivity data in terms 
of power-law scaling^, then the inescapable conclusion is that either the Bell experiment 
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erroneously discarded the data in the true critical region, or the Karlsruhe experiment 
mistook spurious effects for the critical behavior. If this was the case, then a careful scaling 
analysis of both data sets and, if necessary, new and more accurate experiments, should 
be able to show that the conductivity, and possibly other quantities, show scaling behavior 
in one region but not in the other, thus settling the issue. There is, however, another 
possibility. If the theoretical suggest ionl!! that the AMT has random-field aspects is correct, 
then one would expect glass-like features and unconventional scaling similar to what has been 
predictecSl and observedii in classical random-field magnets. 

In this paper we explore these possibilities. In Sec. [TT]we assume conventional scaling, and 
check whether a scaling analysis of the existing data can settle the disagreement between the 
experimentalists. We find that it cannot, mostly due to an insufficient temperature range and 
the lack of precision experiments that probe the critical behavior of more than one quantity 
in a given material. Our analysis suggests, however, various experiments that might be able 
to tell which, if any, of the two doping regions that have been suggested to be the critical one, 
displays conventional scaling behavior. In Sec. [H| we assume instead that the AMT features 
activated scaling of the type found in random-field magnets, appropriately modified for a 
quantum problem. Accordingly, we first develop a general scaling description of a quantum 
glass, and then work out predictions for the critical behavior of various observables. We 
check to what extent the existing data are consistent with these predictions and propose 
new experiments to further investigate this issue. 



II. CONVENTIONAL SCALING 



A. Homogeneity laws 



Let us recall the homogeneity laws for the tunneling or single-particle DOS, N, and the 
conductivity, a, that express conventional power-law scaling,0 

N(t,T) = b~P /v N(tb 1/v ,Tb z ) , (2.1) 



cr(t,T) =b- s/v a(tb 1/u ,Tb z ) . (2.2) 

Here t = n/n c — 1 denotes the dimensionless distance from the critical point, and T is the 
temperature. Since in a quantum problem temperature and frequency scale the same way,0 
one obtains the same homogeneity laws at T = with T replaced by lu, where to is the 
external frequency in the case of the conductivity, and the distance in energy space from 
the Fermi level, i.e. the bias voltage in a tunneling experiment, in the case of the DOS. (3 
and s are the critical exponents for the DOS and the conductivity, respectively, v is the 
correlation length exponent, z is the dynamical critical exponent, and b is an arbitrary length 
scale factor. Analogous homogeneity laws can be written down for all other quantities of 
interest.0 Here we focus on the DOS, since it is the order parameter for the AMT according 
to our recent theoryB~i and since it is easily measured, and on the conductivity since it is the 
most obviously interesting observable in the context of a metal-insulator transition. Note 
that throughout this paper we ignore the possibility of significant corrections to scaling. 
Putting b equal to the correlation length, b = £ ~ t~ u , we obtain 
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N{t,T)=t p F N {T/t vz ) , (2.3) 

and 

a(t,T)=t s F a (T/f z ) , (2.4) 

where F N and F a are scaling functions. 

These homogeneity laws predict that N and a are functions of a particular combination 
of their two arguments t and T, and they have been derived from the traditional description 
of the AMT.0 They also follow from a perturbative treatment of the order parameter theory 
put forward recently.!!! As has been mentioned in these references, it is likely that in the 
second case the perturbative results are misleading, and that the dynamical scaling in the 
physical dimension d = 3 is of activated rather than of power-law type. We will explore the 
consequences of activated scaling in Sec. [TTT] below. For now we assume that the power-law 
scaling expressed by Eqs. ( p73| , |2j4j) is correct asymptotically close to the transition observed 
in Si: P. This could be due to either the OP theory of the AMT not being applicable to Si:P, 
or to the perturbative analogy between the AMT and classical random field magnets being 
misleading. In this section we investigate the experimental consequences of this assumption. 
In doing so it is important to realize that the theoretical values of the critical exponents are 
not known.il 

Another quantity of interest is the correlation function of the local, unaveraged DOS, 

C(t,T;x-y) = (iV(x)iV(y)> , (2.5) 



which according to Refs. is the order parameter susceptibility for the AMT. In Eq. ( |2.5|) , 
iV(x) is the local DOS at the Fermi energy, and < . . . > denotes the impurity average. The 
Fourier transform of C obeys the homogeneity law 

C(t, T; q) = b 2+e -^ C(tb 1/u , TV- qb) . (2.6) 

Here r\ is the usual critical exponent that governs the spatial dependence of the OP sus- 
ceptibility, and 9 > is minus the scale dimension of a dangerous irrelevant variable in the 
random-field problem.!'! In the presence of random-field effects, the perturbative value of 9 
is 2, while in the absence of random-field effects one has 9 = 0. Even though the perturbative 
result cannot be correct, at least not in d = 3, it is very likely that the presence of 9 will 
overcompensate r], and make the scale dimension of C larger than 2. This leads to a strong 
divergence of C(t = OJT = 0;q —>■ 0), the analogue of which has been observed in classical 
random-field magnets.Ba Alternatively, one can consider the homogeneous correlation to find, 

C{t, T;q = 0) = T-^-^l* F c (T/t uz ) , (2.7) 



with F c a scaling function. Equation (|2.7|) predicts a strong divergence as T — > 0, the 
presence of which can be checked experimentally. We will come back to this point. 



B. Scaling analysis of the Bell experiment 

Let us now check whether the data of Ref. [10] are consistent with the above homogeneity 
law, Eq. (|2.4j ). From a double logarithmic plot of the extrapolated T = conductivity 
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the experimentalists determined the critical value of the stress, S c , in their stress-tuning 
experiment to be S c = 6.5 ± 0.2 kbar, and the conductivity exponent s ~ 0.5. From the In- 
dependence of a at a stress value estimated to be close to the critical one, they also inferred 
the value of vz ~ 2.7. In Ref. || it was shown that with these exponent values the data do 
not obey scaling. It must be emphasized, however, that the error in the value of vz is at 
least 3O%.0 Accordingly, let us keep S c = 6.5 kbar and s = 0.5 fixed, but let vz float freely 
to produce the best scaling plot. The result is shown in Fig. [TJ, which replots data from Fig. 
1 of Ref. in a way suggested by Eq. ( p.4j) . With vz = 2.13 one obtains a scaling plot that 
is better than the one in Ref. ^ although its absolute quality is not very good. 

We next check whether the quality of the scaling plot can be improved by changing S c . 
The largest value of S c that is consistent with the error bars given in Ref. |10| is S c = 6.7 kbar. 
Indeed, inspection of Fig. 1 in Ref. [LCj shows that at S = 6.59 kbar the data still show 
the curvature at the lowest temperatures that the authors considered characteristic of the 
insulating regime. Let us therefore assume that the next higher stress value, S = 6.71 kbar, 
was the critical one. With this value for S c , we found that the best scaling plot is achieved 
with s = 0.29 and vz = 1.82, which is shown in Fig. 0. The quality of the scaling plot is 
now much better. 

As Figs. Ill show, the dynamical scaling plot favors a large value of the critical stress, 
S c , and a correspondingly small value of the conductivity exponent s. This requires some 
comments in the light of the determination s = 0.51 ± 0.05 in Ref. [10. The determination 



of S c and s in Refs. p-10 was guided by the desire to achieve a good static scaling plot, 



i.e. a straight line of logo" vs. logi, over as large a t-interval as possible. While this is a 
legitimate and often used criterion for determining the critical point, in the case of Si:P it 
leads to a rather peculiar result: It is found that an exponent of s = 0.5 fits the behavior of 
o~(T = 0) very well out to t « 1.0. While this is remarkable and may well have interesting 
(and presently unknown) reasons, it is very unlikely that the critical region in Si:P is that 
large. Indeed our dynamical scaling plot, which is a more sophisticated test of scaling than 
the static one, shows that it is not. On the other hand, a static scaling plot over a more 
restricted t-range, viz. the data from Fig. 1 of Ref. [T^ (which is the same data set that was 
used to produce the dynamical scaling plots in Figs. ^ and 0), shows that S c = 6.71 kbar and 
s = 0.29 fits the data close to the transition very well, as shown in Fig. |3[. This interpretation 
suggests a size of the critical region of about 1%, which is comparable with the corresponding 
value for most thermal phase transitions.il 

We conclude that the Bell experiment allows for a reasonably good dynamical scaling 
plot, consistent with a conventional power-law scaling interpretation, provided that the 
location of the critical point is adjusted upward, and the value of the conductivity exponent s 
downward, compared to the values given in Ref. [H| The resulting small value of s aggravates 
the problem that results from the inequality v > 2/3 in conjunction with the exponent 
relation s = v. This interpretation is therefore only feasible within theories that allow for 
s 7^ v, as e.g. the order parameter description of Ref. HH Any theory that yields power-law 
scaling with s — v, like those reviewed in Ref. Ill is inconsistent with the Bell experiment 
unless there are large corrections to scaling.!! For later reference we also note that the 
dynamical scaling observed in this experiment is restricted to a small dynamical range of 
just a bit over one decade: For the plot in Fig. ||only data at temperatures T < 60 mK were 
included. Data at higher temperature do not scale, as can be seen from the inset in Fig. ^. 
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C. Scaling analysis of the Karlsruhe experiment 



For the Karlsruhe data a dynamical scaling plot has been given by Stupp et al. in Ref. [T3 . 
For a direct comparison with the Bell data, we have digitized the data from Fig. 1 of Ref. [L2|, 
and plot them in Fig. ^ in the same way as the Bell data in Figs. [I] and |[ We have included 
data for T < 160 mK, but have left out the lowest temperature points on some samples that 
showed obvious rounding effects. We will discuss this rounding in the next subsection. We 
have assumed a critical phosphorus concentration of 3.52 x 10 18 cm -3 , and exponent values 
s = 1.3, and vz = 2.7. This yields the plot shown in Fig. f|. The slight differences between 
our dynamical scaling plot and that of Stupp et al. (who found the optimal value of vz to 
be 3.5) are due to Stupp et al. optimizing over a larger range of t-values, and possibly due 
to some errors introduced by redigitizing the data. Overall, however, the two plots are of 
comparable quality. Fig. [5] shows a static scaling plot analogous to the one shown in Fig. |3|. 
Again, the assumed value of s fits the T = conductivity well over one decade of t. 



D. Discussion of the conventional scaling interpretation 



Sections [II B , f\ C| and the accompanying figures show that the data from both the Bell 



and the Karlsruhe experiment allow for dynamical scaling plots of equal and satisfactory 
quality, even though their results are mutually inconsistent. A check for dynamical scal- 
ing is therefore not sufficient to distinguish between them, and it is necessary to consider 
addditional experimental information to settle the issue. 

At the heart of the discrepancy lie the different values of n c , which differ by 6% between 
the two groups. An independent measurement of n c with an error of about 1% or less would 
therefore suffice to rule out at least one of the two interpretations. Unfortunately, both 
groups claim to have done just that, but again do not agree on the results. Paalanen et al.0 
have measured the dielectric polarizability and in at least one sample have found insulating 
behavior for values of n as little as 1% below their value of n c . Lakner and von LohneysenEd 
have found the thermopower to exhibit a metallic characteristic for values of n above their 
n c , but below that of the Bell group. 

Clearly, what is needed in this situation is an independent determination of n c . The 
existing measurements of the spin susceptibility^ and the specific heat0 are not suitable for 
this purpose since they showed singular behavior in the metallic phase far from the critical 
point. One therefore expects these thermodynamic susceptibilities to show a superposition of 
critical behavior and some noncritical, but nevertheless singular, background, which makes 
them unsuitable for the present purpose. The most obvious observable that should be free 
from such complications, and that is easy to measure, is the DOS. If the conventional scaling 
scenario is correct, then the DOS as a function of t and T should show dynamical scaling, 
as expressed by Eq. (|2.3j), and the quality of the dynamical scaling plot should be equal to 
that of the conductivity with the same parameter values. 

Even such an additional measurement, however, may turn out to be inconclusive unless it 
is carried out at sufficiently low temperatures, and unless the nature of the sample-to-sample 
fluctuations that were observed in the 'tail' region in Refs. HP is clarified. The Karlsruhe 
experiment did not observe such fluctuations, and seemed to yield a smooth conductivity 
as a function of doping. However, this may be misleading since the lowest T reached in 
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this experiment was 62.5 mK, while the sample dependence in the Bell experiment became 
obvious only at lower temperature. In this context it is interesting to note that in the latter 
there is a clear break in the temperature dependence of the conductivity around T = 60 mK, 
see Fig. 1 of Ref. ^ and the inset in our Fig. [|. Furthermore, in the same temperature 
range sample-dependent problems did start to arise in the Karlsruhe experiment, which 
were attributed to thermal decoupling, see Fig. 1 of Ref. [12|. Finally, we note that all of 
the work done by the group at CUNY, both on Si:Piil and on Si:B@, which also yielded a 
smooth behavior of the conductivity, was at temperatures higher than 60 mK. 

We conclude that the existing experimental data on Si:P provide evidence that close 
to the critical point (whose location is only imprecisely known), and at temperatures 
T < 60 mK there are large sample-to-sample fluctuations, and possibly equilibration prob- 
lems. While the data at higher temperatures, and the low-temperature data with the 'tail' 
region discarded, are not inconsistent with conventional scaling behavior, this poses the 
question whether the critical behavior at the AMT might be more exotic than is suggested 
by Eqs. ( |2.1|) - (OP- Indeed, our recent work on an order parameter description of the 
AMT shows that the AMT has random-field aspectsBi Classical random-field magnets are 
well known to displavglassy behavior with exponentially long equilibration times, uncon- 
ventional scaling, etcBa It is natural to assume that similar phenomena can characterize the 
AMT. While any conventional scaling interpretation of the AMT in Si:P will necessarily 
imply that one of the disagreeing experimental groups made a gross error in the determi- 
nation of the critical concentration, an interpretation in terms of random-field physics has 
the potential for explaining the unusual features observed, and the disagreements between 
the experimentalists, in terms of real physical effects that are germane to the AMT. We 
consider this a very appealing possibility. In the remainder of this paper we therefore take 
this suggestion seriously, develop it, and then come back to a discussion of the experimental 
situation. 



III. QUANTUM GLASSY BEHAVIOR, AND ACTIVATED SCALING 

A. Scaling theory of a quantum glass transition 

An important characteristic of a glass transition, as opposed to an ordinary phase tran- 
sition, is the occurrence of extremely long time scales. While critical slowing down at an 
ordinary transition means that the critical time scale r grows like a power of the correlation 
length, r ~ £ z , with z the dynamical critical exponent,^ at a glass transition the critical 
time scale grows exponentially with £,@ 

In (t/to)~^ , (3.1) 

with t a microscopic time scale, and if> a generalized dynamical exponent. As a result of 
such extreme slowing down, the system's equilibrium behavior near the transition becomes 
inaccessible for all practical purposes: Even at the smallest feasible frequencies (i.e., the 
inverse of the longest feasible waiting times, say, days), finite frequency effects become 
noticable ("the systems falls out of equilibrium") at modest values of £ or t. It has been 
proposed that the phase transition in classical random-field magnets is of this type, and 
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experiments have confirmed this conjecture.il The physical picture behind this model of 
random-field magnets is as follows. The frustration induced by the competition between the 
exchange interaction and the random field leads to large clusters of misaligned spins, that 
is, locally ordered spins that 'point the wrong way', within the ordered phase. Even though 
aligning these clusters leads overall to a lower free energy, it requires a large free energy 
barrier to be overcome. These free energy barriers grow like as a function of some length 
scale L, and near the critical point they diverge like the correlation length £ to the power ip. 
The exponent ifi is therefore often referred to as the 'barrier exponent'. Via the Arrhenius 
law, this leads to Eq. 

In a quantum system one expects time and inverse temperature to show the same scaling 
behavior, irrespective of whether the critical slowing down follows an ordinary power law, 
or Eq. ( p.l| )B Quantum mechanics thus makes it even harder to observe the static scaling 
behavior, since in addition to exponentially long times or small frequencies it requires ex- 
ponentially small temperatures as well. Under realistic experimental conditions the system 
will either fall out of equilibrium, or pick up finite temperature effects. This is a crucial 
point which we will have to keep in mind for any discussion of experimental consequences 
of our ideas. 

The role played by temperature in a glassy quantum system can be seen explicitly in 
Fisher's recent study of a quantum-mechanical Ising spin chain in a transverse random 
magnetic field.@ This system is closely related to the classical McCoy- Wu model, for which 
a number of exact results have been obtained. EH Since the AMT has been shown to be 
a quantum phase transition with random-field aspectsi we believe that these results are 
qualitatively relevant for our purposes, and will often use them for comparisons. The physical 
idea analogous to the one explained above for spin systems is that while a repulsive electron- 
electron interaction always leads to a decrease in the local DOS, the random potential can 
in general lead to an increase in the local DOS as well. The competition between these two 
effects leads to frustration and to, for example, large insulating clusters within the metallic 
phase. Delocalizing these large clusters requires energy barriers to be overcome which again 
grow like ^ as the transition is approached. 

Another question that arises is the scaling behavior of an external magnetic field. In 
a spin system the magnetic field, just like the temperature, obviously sets an energy scale 
and thus should depend exponentially on the length scale. In an itinerant electron system, 
on the other hand, the magnetic field plays a dual role: It influences the orbital motion 
of the electrons, and in that capacity it acts like a length. However, it also couples to the 
electron spin via the Zeeman term, and there it acts like an energy. We therefore expect the 
slowest dependence of a given observable on the magnetic field to be a logarithmic one, with 
power-law corrections. A further consequence of the barrier model is that the frequency or 
temperature argument of scaling functions is expected to be ln(r/r )/ ln(T /T), rather than 
tT as in Eqs. ( |2.1| , |2.2| ). The reason is that one expects a very broad distribution of energy 
barriers and hence of relaxation times r. The natural variable is therefore In r rather than 
r.0 This is often referred to as activated scaling. 

We next face the question of which observables can be expected to obey homogeneity 
laws analogous to the ones given by Eqs. ( |2.1| , |2.2| ). First of all, we have to remember 



that the homogeneity laws hold for averaged quantities, with the average including both a 
quantum mechanical and an impurity or ensemble average, and that we have to distinguish 
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between self-averaging and non-self averaging quantities.il For the former, their probability 
distribution for the ensemble average becomes normal with a vanishing variance in the 
thermodynamic limit, so their average value coincides with the most probable or typical 
one. For the latter, either the probability distribution remains broad (e.g. log-normal), 
or the most probable value is different from the average one, or both. If the observable in 
question is non-self averaging due to a log-normal distribution, then one expects its logarithm 
to be a self-averaging quantity. 

It is well known that in a system with quenched disorder the free energy is self-averaging, 
while the partition function is not, and correlation functions in general are not, either Jl§Ea 
Therefore, all thermodynamic quantities, which can be obtained as partial derivatives of 
the free energy, are self-averaging. For a general thermodynamic quantity, Q, one might 
therefore expect a homogeneity law 

W,T)=r-F Q _|_) , (3.2, 

where xq is the scale dimension of Q, Fq is a scaling function, and T is a microscopic 
temperature scale, e.g. the Fermi temperature. A finite frequency will have the same effect 
as a finite temperature; as mentioned above, frequency and temperature are interchangeable 
in a scaling sense. A magnetic field dependence will be added later, when we discuss the 
magnetization. 

It is obvious, however, that Eq. ( [3.20 can not hold for all self-averaging quantities, or even 
for all thermodynamic ones. Suppose Q is some self-averaging quantity, and suppose Eq. 
( [3.20 holds for Q. Then P = Q/T is also self- averaging, but since Fq is not a homogeneous 



function of T, Eq. ( |3.2|) will not hold for P. While this is a rather trivial 'breakdown of 
scaling', it has observable consequences as we will see below. Here we assume that whether 
or not a given thermodynamic quantity obeys Eq. ( |3.2| ) can be decided by dimensional 



analysis: If the quantity contains scale dimensions of time or energy, then it does not obey 
Eq. ( p. 20 , else it does. For instance, the free energy density can not be written in this form, 
while the entropy density can, etc. However, the question whether or not a particular energy 
acts like an inverse time in a scaling sense can be a nontrivial one. We will further consider 



this point when we explicitly discuss various observables in Sec. |III B| below. 

A rather different class of observables is formed by transport coefficients, like, e.g., the 
charge or heat diffusivity. Since they are directly related to a relaxation time, we expect 
them to be non-self averaging, while their logarithms should be self-averaging. Let S be an 
unaveraged tranport coefficient, i.e. its value for a particular sample or impurity arrange- 
ment. Then we expect its logarithm to obey 

.„E)( ( ,T) = ^(«V % _|_) , (3.3) 



with Fs a scaling function. Notice that Eq. ( |3.3|) describes only the leading, i.e. logarithmic, 
scaling behavior, and neglects power-law corrections to scaling. The 'scale dimension' of 
(\n Ej is necessarily plus or minus ip, with the sign depending on whether the quantity 
vanishes or diverges at the transition. 



As with conventional scaling at ordinary phase transitions, Eqs. (BT2L |3.3| ) hold only for 



the singular parts of the respective quantities, and in general there will be nonvanishing, 
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analytic background contributions. In the case of Eq. ( |3.3| ) another complication is to be 
expected. A general transport coefficient S is related, by means of an Einstein relation, to 
the corresponding diffusivity A via H = xA, with \ an unaveraged susceptibility. Since x 
is a thermodynamic quantity, ln% is not expected to show scaling behavior. If \ is critical, 
one therefore expects a critical, non-scaling background contribution to AnH^, in addition 
to the scaling part given by Eq. ( |3.3| ). We will come back to this. 

In the above paragraphs we have stated all of the assumptions that enter our scaling 
theory of a quantum glass. In the remainder of this paper we explicitly discuss the behavior 
of a number of specific observables that are of interest in the context of the AMT. We work 
out the consequences of our assumptions, and compare the results with those obtained from 
conventional scaling, Sec. O, and with the experiments on Si:P discussed above. 

B. Discussion of observables 

In this section we discuss explicitly the consequences of our scaling assumptions for 
various observables. We start with thermodynamic quantities, for which Eq. (|3.2|) and 
the related discussion above are relevant. Then we turn to the electrical and the thermal 
conductivity as examples of transport coefficients which realize the scaling behavior shown 
in Eq. Q. 

1. Tunneling density of states 

Let us first discuss the tunneling density of states, N, as a function of t and T. We 
restrict ourselves to the density of states at the Fermi level, since according to the discussion 
in connection with Eq. ( |3.2|) a finite frequency or energy, which measures the bias voltage or 
the distance from the chemical potential, will have the same effect as T. In the theory put 
forward in Refs. |]-|6], N is the order parameter for the AMT, and its scale dimension follows 
directly from the order parameter field theory to be xjq = d — 9 — 2 -\- rj. Comparison with 
Eq. (|2.1|) gives the exponent relation!! 

P=^d-9-2 + r] ) . (3.4a) 

From Eq. ( |3.2| ) we then obtain the generalization of Eq. ( |2.1D to the case of activated scaling 
as 

N(t, T) = b-V"F N (tb^ , ln( jf /T) ) . (3.4b) 

Here is a scaling function, and 9 is an exponent related to a dangerous irrelevant vari- 
able that is characteristic of the random-field problem.^ In a classical random-field problem 
thermal fluctuation are dangerously irrelevant, and 9 is minus the scale dimension of temper- 
ature. In the present context 9 expresses the fact that quantum fluctuations are dangerously 
irrelevant at the AMT.i 

A fundamental question that arises in this context is how many independent exponents 
are needed to describe the AMT. In a classical random-field problem there are three inde- 
pendent static exponents, e.g. u, rj, and In the AMT theory developed in Ref. |5]|| it 
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turned out that the dynamical exponent z was not independent, reflecting the irrelevancy 
of quantum fluctuations, so that there were still three independent exponents. Here if> has 
taken over the role of z, and the question is whether or not it is independent. In order 
to decide this, let us recall the physical meaning of the two exponents 9 and if). As noted 
above, 9 is related to a dangerous irrelevant variable (DIV), u, which vanishes as a func- 
tion of length scale L like u ~ L~ e . The random-field fixed point is characterized by uA 
scaling to a constant, with A the random potential energy scale.0 Hence A must diverge 
as L e . The free energy landscape of the random-field problem is a complicated one, with 
many near-degenerate valleys that are separated by energy barriers with saddle points. As a 
function of length scale L, one expects a typical valley elevation to be related to the random 
potential and therefore grow like L e . The typical saddle point elevation grows like L^, which 
defines the exponent if). In order for this picture to be consistent, one must have if) > 9, as 
pointed out by Fisher in Ref. |38|. Classically, only the valleys contribute to the free energy, 
so one expects the exponent characterizing the dangerously irrelevant thermal fluctuations 
to be 9, and the barrier exponent to be if), and in general the two will be independent. 
Quantum mechanically, however, the saddle points also contribute to the free energy, and 
one cannot distinguish between the barrier exponent and the exponent that expresses the 
fact that quantum fluctuations are dangerously irrelevant. We therefore expect 9 — if> in 
the quantum case, although in our notation we will continue to distinguish between the two 
exponents. This leaves us with three independent exponents, e.g. v, rj, and if> — 9. A fourth 
one will be necessary when we discuss sytems in external magnetic fields in Sees. |l [ B 4| and 



[jrBJ below. 

As in the case of ordinary power-law scaling, we can eliminate the arbitrary parameter 
b from Eq. ( |3.4b| ) and write N(t, T) in two different scaling forms to emphasize either the 
static or the dynamic aspects of the scaling law. The crossover between the two types of 
scaling occurs at a temperature T x for which the two arguments of Fn are equal. Ordinarily, 
this criterion would lead to a power-law dependence of T x on t, but activated scaling implies 
T x ~ exp(— l/t u ^). As a result, the static scalingregion will be very small unless uif) is ver 
small. Since v is bounded from below, v > 2/rf,B3 this would require if> to be very small.c 
Whether or not the static scaling behavior, N(t, T) « N(t, T = 0) ~ t" is observable will 
then strongly depend on the precise value of vif), on the size of the critical region, and, 
to a lesser extent, on the value of the microscopic temperature scale T . We thus put 

= ln(T /T), and write Eq. ( ggg ) as 

N(t, T) = - - G N [f* \n(T /T)] , (3.5) 



3 



MVO] 



The scaling function is related to the function F N in Eq. (|3.4b| ) by Gat(x) = F^{x l ' v ^ , 1), 



and has the properties G^{x — > oo) ~ x®l v ^ , and Gn(x — >• 0) — > const. 



Equation ( |3.5| ) makes a qualitative prediction that can be used to check experimentally 
for glassy aspects of the AMT: Measurements of the tunneling density of states very close 
to the transition should show an anomalously slow temperature depedence, i.e. iV should 
vanish as some power of InT rather than as a power of T. While this is a straightforward 
check in principle, in practice it may require a very large T-range to distinguish between the 
two possibilities. For instance, in classical magnets the frequency had to be varied over seven 
decades in order to convincingly demonstrate the presence of activated scaling.il However, 
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measurements over a smaller dynamic range would also be of interest, since they would 
put experimental bounds on possible values of ip. This is particularly important since in 
the absence of any information about the value of if) it is impossible to tell whether at a 
given temperature one is in the static or the dynamic scaling regime. Unfortunately, to our 
knowledge all measurements of N close to metal-insulator transitions have been performed at 
fixed (and rather high) temperatures, so that no T-dependent data are available for analysis. 



2. Order-parameter susceptibility 

We now turn to fluctuations of the order parameter. We first give a statistical argument 
that N is indeed, as assumed above, a self-averaging quantity. This also sheds some light 
on the crucial role played by the electron-electron interaction in our theory for the AMT. 

Let us write the unaveraged order parameter, N, as its average plus fluctuations, N = 
N + 5N, and consider the mean-square fluctuation ((5iV) 2 ). At T = in a system of size 
L one has 

(5N) 2 ) = N$(t,L) , (3.6) 
with $ some function of t and L. Far away from the critical point we expect $ = const, 

1 /2 

which leads to as usual. In the critical region, on the other hand, 

we expect $ to scale, 

${t,L) = b' y/v <$>{tb 1/v ,Lb- 1 ) , (3.7) 

with 7 the critical exponent for the order parameter susceptibility. At criticality, Eq. 
implies (|(<5iV) 2 ^ ~ N L 2 ~ v ~ L d+2- '', where we have used the exponent relation 7 = v(2—rj), 
and N ~ L d . For the root mean-square order parameter fluctuations this means 

(((5N) 2 ) /N 2 ) 1 ' 2 ~ L^ d -^ 2 ~ L~ 2P ' V . (3.8) 



For the last relation in Eq. ( |3.8| ) we have assumed hyperscaling to be valid. In its absence the 
argument needs a trivial modification, but the end result is still given by the far right-hand 
side of Eq. (pp. 

Equation ( |3.8|) implies that order parameter fluctuations at the critical point become 
small in large systems provided 7] > 2 — d, or /3 > (the first condition depends on hy- 
perscaling, while the second one does not). This has some interesting consequences. In the 
noninteracting localization problem, one has 77 = 2 — d, and (3 = O.cffl As can be seen from 
the above discussion, this means more than simply that the order parameter is uncritical 
in the localization problem: It indicates that the fluctuations of the density of states are 
independent of the system size, and as large as the average. Consequently, one expects that 
the density of states in a system of noninteracting disordered electrons has a very broad 
distribution, and that the Anderson transition is pathological from a Statistical Mechanics 
point of view. All of this is consistent with explicit studies of the Anderson transition.il Our 
order parameter description of the AMT, on the other hand, leads to (3 > 0. Therefore the 



density of states will be self-averaging in accord with our assumptions in Sec. [Ill A| above, 
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and there are no obvious obstacles for the description of the AMT in terms of the standard 
concepts for continous phase transitions. 

We now discuss the order parameter correlation function, which we define as 

C(t, T; x - y, lu) = (iV(x, e F + u/2) N(y, e F - to/2)) , (3.9) 



or its spatial Fourier transform, C(t,T;q,u). Here iV(x, e) denotes the unaveraged, local 
density of states at energy e. From an analogy with Eq. Q2.6D one expects C to scale, and 
to show an anomalously strong divergence as q — > 0. However, in random- field systems 
the correlation function at nonzero wavenumber is likely to not be normally distributed,^ 
while at q = no such complications are expected to occur. Moreover, in our quantum 
system a small-g divergence would be cut off by a finite temperature. We therefore restrict 
ourselves to a discussion of the homogeneous correlation function, i.e. the order parameter 
susceptibility, which obeys 

C(«,T i8 = 0, W ) = ^F (* . (3 . 10 ) 

At criticality, the order parameter susceptibility diverges as the temperature goes to zero, 
but only as a power of In T, 

C(t = 0,T;q = 0,oj = 0)~[ln(T Q /T)f +e ~ r > W . (3.11) 

The local density of states is measurable (on a perfect surface) with a scanning tunneling 
microscope, or STM. By measuring the density of states both by means of a tunnel junction 
and by means of an STM one can therefore check whether it is indeed a self-averaging 
quantity: If it is, then in a large system fluctuations of the local density of states should 
be small, and both the local and the junction measurements should give the same result. 
Furthermore, by measuring N(x) across a sample at pairs of points with a fixed separation 
it should be possible to measure the correlation function C, and to check the prediction of 

Eq. (fmD. 



We conclude this subsection by noting that C, at least away from the AMT, contains 
an uncritical, 'mesoscopic' power-law singularity due to hydrodynamic diffusion modes.0 In 
order to distinguish this singularity from the one given by Eq. (|3.11| ) one can use arguments 
like those used for the classical random field problem with a conserved order parameter. El The 
key idea is to consider a small but finite wavenumber, q, subject to a number of constraints. 
First we require a self-averaging quantity which according to the discussion above Eq. ( |3.10| ) 
requires g(< 1. Second, we require the critical part of C to effectively be at q = which 
leads to the restriction £q <C 1/ (\nT /T) 1 ^ , with £ a microscopic length on the order of 
the inverse Fermi wavenumber. Third, the noncritical mesoscopic contribution should be T- 
independent. This will be the case if q£ ^> (T/T Q ) 1//2 , since this singularity is due to diffusion. 
It can be readily verified that these three conditions can be simultaneously satisfied. 



3. Specific heat 

Let us now consider the entropy density, s(t,T) = df/dT, where / is the free energy 
density. As a thermodynamic quantity, s is expected to be self-averaging, and since its 
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dimension is that of an inverse volume (in units chosen such that ks = 1) we can write for 
the singular part of s, 

s i t,T) = m[t^,^ m ) . (3.12) 

The specific heat is obtained from s by means of a logarithmic derivative with respect to T, 
cy = ds/dlnT, which yields for the singular part of cy, 



c v (t,T) =b- {d - e+ ^ F cv [tb 1/u . 



ln(T /T) 



[HTo/n 



l+(d-6)/il> 



G CV ((T/T f") . (3.13) 



Here the scaling functions G Cy and F Cv are related by G Cv (x) = F Cv ((— lnx) 1 /^, l). 

The fact that the specific heat must vanish at T = puts constraints on the possible 
behaviors of the function G Cv at small values of its argument. One possibility is G Cv (x — > 
0) — > const. However, a more natural assumption is that G Cv (x — > 0) vanishes like a power 
of its argument. This possibility is realized, e.g., in the model studied recently by Fisher.il 
In that case, the specific heat in the critical region goes like 



const Xf* 



C " (i ' r)= [ln(rI/T)]^-»>/* ' (3 - M) 



If we substitute the exponent values that are appropriate for the model of Ref. viz. 
u = 2, ip = 1/2, 8 = 0, d = 1 and const = 2, then we recover from Eq. ( |3.14 ) Fisher's result 
for the transverse random-field Ising chain. The most interesting aspect of this result is the 
continuously varying exponent in the numerator, which leads to successive derivatives of 
cy becoming divergent as the critical point is approached. This behavior means that there 
is a Griffiths phase, or rather Griffiths region, away from the critical point within which 
certain observables become divergent at various values of t.0 Regardless of the behavior of 
G Cv (x — > 0), we have a non- Fermi liquid behavior of the system in a finite region around 
the critical point: The specific heat coefficient j(t,T) = Cv(t,T)/T diverges as T — > 0, even 
away from criticality. 

For a discussion of the experimental consequences of Eqs. ( 3.13 , 3.14 ) one must keep 
in mind that the singular contribution to cy is additive to the noncritical Fermi liquid 
background that is linear in T, and that the singular part will dominate only at sufficiently 
low temperatures. Measurements of the specific heat in SiiPEl'Ill have indeed observed non- 
Fermi liquid behavior both near the critical point, and rather far away from it in either phase. 
This has been interpreted in terms of local magnetic moments.El The relation between the 
local moment and quantum glass pictures is currently unclear. It is interesting to note that 
both predict singular behavior of thermodynamic quantities away from the critical point, and 
it is conceivable that the glass picture is related to an interacting local moment description. 
Also, from the discussion in Sec. |I| we suspect that the lowest temperatures reached in the 
Si:P experiments (~ 30 mK) were not low enough to be in the critical region. Since going 
to substantially lower temperatures is not realistic, it would be desirable to have similar 
measurements performed on a system with a higher Fermi temperature than Si:P, which has 
Tp w 100K near the critical P concentration. 
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4- Magnetization 



In order to discuss the magnetization, and the magnetic susceptibility in the next sub- 
section, we need to add an external magnetic field H to our discussion. As discussed in Sec. 



Ill A| , the leading effect of a magnetic field will come from its coupling to the electron spin, 
and will scale the same way as the temperature does, viz. T ~ exp(fe^), H ~ exp(6^), both 
up to multiplicative power-law corrections. In general one therefore expects T/H ~ 6*^. 
with an exponent that characterizes differences in the corrections to scaling of T and 
Dimensionally, the magnetization m = df/dH is an inverse volume times a temperature 
divided by a magnetic field, and therefore the scale dimension of m is d — 9 — (pip. Therefore 
we have, 

m(t, T, H) = b- d+e+ ^ F m tb 1/u , — — — , (H/T)b^ , — — -A -—- . (3.15) 

v ' ' > m \ ' ln(T /T) V 1 ' ' ln(T /(T+ Hb^)) J V ; 

Here the last argument of the scaling function F m expresses the fact that the Zeeman energy 
provided by the magnetic field helps the system to overcome free energy barriers, and its 
functional form is motivated by the fact that the effect of a nonzero H will always be cut 
off by a nonzero T. However, the reverse is not true, which is why one still needs the second 
argument containing only T. The third argument contains the physics due to fluctuations 
within a given free energy valley, with no attempts to climb over barriers. 

The most interesting consequence of Eq. ( |3.15|) is the leading if-dependence of m at 
T = at criticality, which is 

m (^0,r = 0,g)~ Hvg) |_^_ w , (3.16) 

where we have assumed that F m (0, 0, oo, 1) is a finite number. If we substitute the exponent 
values that are appropriate for the model of Ref. [35|, namely cp = (1 + Vo)/2 and d, 6, and 



ip as given after Eq. ( |3.14| ), then we recover Fisher's result for the transverse Ising chain. 



The physical interpretation of Eq. ( |3.16|) is as follows. Equation (|3.15|) says that the 



magnetic degrees of freedom are glassy, and relax slowly just as the singlet or DOS degrees 
of freedom. If the system is cooled in a magnetic field then as the magnetic field is turned 
off, the magnetization vanishes slowly as a function of both field and time. This behavior is 
reminiscent of that found in random field magnets. The difference is that here no long range 
magnetic order develops across the transition. However, we do expect that experiments that 
examine the difference between field cooling and zero-field cooling will be very interesting, 
as they are in random-field magnetic systems.^ 



5. Magnetic susceptibility 

So far all quantities we have discussed have been both self-averaging (exept for the 
order parameter correlation function at nonzero wavenumber) and scaling, i.e. they obey 
homogeneity laws of the type given in Eq. ( |3.2| ). The magnetic susceptibility, Xm — dm/dH, 
is clearly self-averaging, but does not obey a homogeneity law for the reasons explained after 
Eq. (|3.2|) . Nevertheless, we can obtain the functional form of \ m by differentiating Eq. (|3.15 ) 
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with respect to H. The leading behavior of the zero-field susceptibility, which is produced 
by the third argument of the scaling function F m , is 

Xm (t,T) = [ln(T /T)n^+ 2 ^ I G Xm ((T/T Q )*+) . (3.17) 

In an ordinary, power-law scaling scenario one would expect the scaling function G Xm to 
behave such that at T = the susceptibility is finite for t ^ 0. Here this is not possible since 
Xm is not a homogeneous function of T. Rather, we conclude that the magnetic susceptibility 
will diverge as T — > in a region of finite size around the critical point. This divergence is 
power-law with logarithmic corrections, and the exponent of the power-law is a continous 
function of t, 

-1 + const xC* 

Xm{t,T) = [HT /T)f- e ^ • (3,18) 

Clearly, this is another manifestation of the Griffiths phenomenon discussed above in Sees. 
|111 B 3| and |111B4| , and again our result is consistent with that of Ref. for the transverse 



random-field Ising chain. 

As mentioned below Eq. ( |3.14 ), there is some similarity between our results and those 



obtained from the local mome"ictureS Both lead to a divergent magnetic susceptibility 
and specific heat coefficient in the metallic phase, in qualitative agreement with experiments 
on Si:P and Si:P,B.llHl§ The main difference is that we predict a critical singularity for either 
quantity, while the local moment picture yields thermodynamic anomalies that are decoupled 
from the AMT. Several points should be kept in mind, however. First, the coupling, or 
absence of it, of local moments to the AMT is an unsolved problem.0 Second, Bhatt and 
FishercO have pointed out that interactions between local moments may considerably weaken 



the effects found in Ref. IS 



6. Density susceptibility 

The thermodynamic density susceptibility, dn/dfi, is not directly measurable in a three- 
dimensional system. However, it is of interest since it enters the Einstein relation between 
the electrical conductivity and the mass diffusion coefficient and can therefore influence the 
critical behavior of the conductivity. 

As a thermodynamic quantity, dn/dfi = <9 2 //<9/z 2 is self-averaging, but it is not a scaling 
quantity in the sense of Eq. ( |3.2|) since / is not. dn/d/j, thus belongs in the same category as 
the magnetic susceptibility, namely that of self-averaging, non-scaling observables. In order 
to determine the critical behavior of dn/d/j,, we first note that the chemical potential \i in 
this derivative does not scale like an energy, but rather like the correlation length to some 
power. This can be seen from the explicit formulation of the order parameter field theory 
for the AMT,! whose only dependence on the chemical potential is in the /i-dependence of 
t. We thus write 

dn/dfi = d 2 f/dfi 2 = (d 2 f/dt 2 ) (dt/dfx) 2 + (df/dt) (dH/dfj 2 ) . 
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At an ordinary quantum phase transition one would expect / ~ &~( d + z ). At glassy quantum 
phase transitions, z effectively diverges so that the free energy does not satisfy a simple ho- 
mogeneneity law. Instead, one expects, schematically, b~ z ~ T ~ exp(— b^) ~ exp(— l/t v ^), 
i.e., the singular part of the free energy has an essential singularity in t at zero tempera- 
ture. This argument is consistent with the results in Sees. |I IB 3| - [III B 5| where we showed 
that the thermodynamics of our model is qualitatively the same as that of the transverse 
random-field Ising chain.0 At T = 0, the latter in turn is equivalent to the model consid- 
ered by Shankar and Murthy,0 who found the singular part of the free energy to behave 
like f(t,T = 0) ~ t 1 /*. We therefore expect this behavior to qualitatively hold in our case 
as well, which means that df/dt and d 2 f/dt 2 vanish exponentially as t — > 0. dt/dfi and 
d 2 t/dfi 2 , on the other hand, can diverge at most like a power of 1/t. This can be seen as 
follows. Let /i c be the critical value of /x at a given value of the disorder, t must vanish as 
/j — > fi c , and it can do so either as a power or as an exponential function of of /i — \x c . In the 
latter case all of the derivatives of t with respect to \x also vanish exponentially, while in the 
former case they may at most diverge like a power of — /i c ). Consequently, the singular 
part of dn/dn at zero temperature must vanish exponentially as t — > 0. This is in contrast 
to the conventional scaling scenario, which yields a power-law dependence of dn/dfi on t in 
the framework of the order parameter field theoryQ and an uncritical dn/dfx in the 2 + e 
expansion. 00 

In addition to this singular part, one expects in general a nonvanishing analytic back- 
ground contribution to dn/dfiE The electrical conductivity will then have the same critical 
behavior as the charge or mass diffusivity. However, if in a particular system, or for par- 
ticular parameter values, that background contribution should vanish, then the resulting 
exponential vanishing of dn/dfi will lead to complications in the scaling description of the 
conductivity, as was mentioned in the discussion of Eq. (|3.3|) above. 



7. Electrical conductivity 

We now turn to the behavior of the electrical conductivity. As explained in Sec. [Ill A| , 
this is not a self-averaging quantity. We consider instead the unaveraged conducitivity, 
a, and define l a =< ln(er /cr) >, with cr a suitable conductivity scale, e.g. the solution 
of the Boltzmann equation. According to Eq. (|3.3|) and Sec. [Ill B 6| we expect l a to be 
self-averaging, and to obey 



l a (t 1 T) = b^F a (tb 1 '\ 



b"' 



HTa/T)J 

= ln(T Q /T) G a (rtn(T /T)) . (3.19) 

This holds if the density susceptibility has a nonvanishing uncritical background contribu- 
tion, as one usually expects to be the case. If dn/d^i vanishes at criticality, then there will 
be a critical, non-scaling background contribution to l a , as explained in connection with Eq. 

(PD. 

As in the case of Eq. ( |3.13| ) or ( 3.17| ), the behavior of the scaling function G a for large 



values of its argument is a priori unclear. However, we can use physical arguments to 
determine it. Let us define S = cr ex P( — la) as a measure of the conductivity. If G a (x — > oo) 
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vanished faster than 1/x, than E would approach <To even for arbitrarily small t ^ as T — > 0. 
Since <Jq is a noncritical quantity, this is unphysical. On the other hand, if G a (x — > oo) 
vanished more slowly than then E(T = 0) would vanish even for t ^ 0. However, for 
t 7^ there are no infinite free energy barriers, and hence density fluctuations are able to 
relax and X must be nonzero. We therefore conclude that G a (x — > oo) ~ 1/x. This yields 

T = 0) ~ exp(-l/t^) , (3.20a) 

and 

E(t = 0,T) ~T Gct(0) . (3.20b) 

Note that at zero temperature E vanishes exponentially with t, and that at the critical point 
E vanishes like a nonuniversal power of T. 

The conclusion that is most important with respect to the interpretation of experimental 
results is that the conductivity, a, is not self-averaging, while In a is self-averaging and 
< In a > scales. While in Sec. we have shown that the existing data for the conductivity 
do allow for scaling plots, we have also seen that many 'strange' features in the experimental 
results, especially in the ultra- low temperature results of the Bell experiment, must be 
ignored in order to reach that conclusion. Also, it was necessary to let all exponents and 
the position of the critical point float. Within our current activated scaling scenario there is 
no reason to believe that a would scale if one obtained bounds on the exponents and on n c 
by measuring thermodynamic quantities (which do scale even under the current scenario) 
at the same low temperatures as the conductivity. < In a > does scale, but would be 
hard to measure. In other words, if activated scaling is present but conventional scaling 
is used for analyzing experiments, then better experiments will make things worse rather 
than better. Furthermore, a is predicted to not be a self-averaging quantity, but to have a 
broad probability distribution. Measurements of a at sufficiently low temperatures should 
therefore show large sample-to-sample fluctuations. 

We propose that the unusual features that were observed in the experiments on Si:P, 
particularly in the ultra-low temperature Bell experiment, and which we have reviewed 
in Sees. | and P], are manifestations of the 'glassy' behavior that we have derived above. 
The fact that the observed anomalies became stronger at lower temperatures is certainly 
consistent with this. The fact that other experimentsB did not provide any indications for 
a not being a well-behaved quantity is not a valid counterargument, since they all stayed 
above, or barely got below, the 60 mK where there is a clear break in the T-dependence of 
the conductivity, see Ref. [II] and Fig. |[ In order to further check this proposal, one should 
measure thermodynamic quantities, preferably the tunneling density of states, together with 
the conductivity at as low temperatures, and over as wide a temperature range, as possible. 
A system with a higher Fermi temperature than doped Si would be advantageous, since it 
would alleviate the need for ultralow temperatures. Ni(S,Se) 2 may be promising in this 
respect. In a very recent interesting experiment, Jin et al. have found that, although 
conductivity data down to T = 30 mK do allow for a conventional dynamical scaling plot, 
there are hysteresis effects which may be indicative of a glass-like behavior of the electrons^ 
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8. Thermal conductivity 

We finally consider the electronic contribution to the thermal conductivity, which is the 
product of the specific heat and the heat diffusivity. For the same reasons as in the case of 
the electrical conductivity we expect the thermal conductivity, k, not to be a self-averaging 
quantity. We define l K = (ln(K / K )), with n the Boltzmann value, and expect0 

« ( ,T)=6*F„( ( ^_|y . (3. 21 ) 

This equation can be discussed analogously to Eq. ( |3.19| ) for the conductivity. If we define 
K = exp(— l K ) as a measure of the thermal conductivity, then we obtain an interesting 
prediction for the generalized Wiedemann-Franz ratio K/H at criticality: 

K(t = 0,T)/£(t = 0,T) ~ T GK ~ Ga , (3.22) 

with G K and G a nonuniversal numbers (see Eq. ( 3.20b|) ). This is in sharp contrast to the 
conventional scaling description of the AMT, which predicts that the Wiedemann-Franz law 
k/o =< k > I < a > ~ T holds even at the transition. @ 



IV. SUMMARY 

We conclude by briefly summarizing the results of this paper. We have employed both 
conventional and activated scaling scenarios to analyze experiments on the metal-insulator 
transition in doped semiconductors, most notably Si:P. Our main goals were to understand 
the discrepancies between different experimental findings, and to work out and analyze the 
glassy dynamical features of the transition that are suggested by recent theoretical advances. 

In Section || conventional scaling ideas were used to interpret existing experimental 
data. The most important conclusions were that existing experiments for the conductivity 
are inconsistent with each other and that, at least in Si:P at very low temperatures, there are 
large sample-to-sample fluctuations, and possibly equilibration problems, sufficiently close 
to the critical point. 

In Sec. |TJ we assumed that the AMT is a quantum glass transition, and we developed a 
general description of such a transition. Our chief results are as follows: (1) The specific heat 
and spin susceptibility are singular as T — > even in the metallic phase, see Eqs. ( p,14| , |3.18|) . 
These results are consistent with existing experiments, and the theory given here provides 
an alternative to the previous explanation in terms of noninteracting local moments. (2) The 
DOS is the order parameter for the quantum glass transition and it is both self-averaging 
and a scaling quantity, see Sec. |I 1 B 1| . At criticality, it is predicted to vanish logarithmically 



with temperature, see Eq. ( |3.5| ). The critical behavior of the OP susceptibility has also been 
discussed. (3) The electrical conductivity, a, is so broadly distributed that it is not a self- 



averaging quantity, but In a is both self-averaging and a scaling quantity, see Sec. |I[_B7. 



This result was used to explain the sample-to-sample fluctuations in a that were observed in 



Si:P. In Sec. [HI] we also suggested a number of additional experiments to test the hypothesis 



that the AMT is a quantum glass transition. 
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FIGURES 

FIG. 1. Dynamical scaling plot of the conductivity data from Fig. 1 of Ref. The plot 
assumes a critical stress S c = 6.5 kbar, and exponent values s = 0.5, vz = 2.13. Only data in 
the temperature range T < 60 mK have been included in the plot, and different symbols denote 
different stress values, from 5" = 6.59 kbar to S = 8.03 kbar. We have chosen To = 100 K, and the 
relation between t and S — S c was taken from Ref. |, viz. t = (S - S c ) 5.4 x 10~ 3 (kbar)- 1 . 

FIG. 2. Same as Fig. |l], but with S c = 6.71 kbar, s = 0.29, and vz = 1.82. The inset shows that 
the data cease to scale once the temperature region 60 mK < T < 225 mK is taken into account. 

FIG. 3. Static scaling plot of the conductivity at T = with the same parameter values as for 
the dynamical scaling plot in Fig. [2|. The line corresponds to an exponent s = 0.29. 

FIG. 4. Dynamical scaling plot of the conductivity data from Fig. 1 of Ref. [l^. The plot 
assumes a critical P density n c = 3.52 x 10 18 kbar, and exponent values s = 1.3, vz = 2.7. Only 
data in the temperature range T < 160 mK have been included in the plot, and different symbols 
denote different P densities, from n = 3.55 x 10 18 cm -3 to n = 3.69 x 10 18 cm -3 . We have chosen 
T = 100 K. 

FIG. 5. Static scaling plot of the conductivity at T = with the same parameter values as for 
the dynamical scaling plot in Fig. |4|. The line corresponds to an exponent s = 1.3. 
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